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ABSTRACT 

Context. Studying the trajectories of objects like stars, globular clusters, or satellite galaxies in the Milky Way allows the dark matter 
halo to be traced but requires reliable models of its gravitational potential. 

Aims. Realistic, yet simple and fully analytical, models have already been presented in the past. However, improved, as well as new, 
observational constraints have become available in the meantime, calling for a recalibration of the respective model parameters. 
Methods. Three widely used model potentials are revisited. By a simultaneous least-squares fit to the observed rotation curve, in-plane 
proper motion of Sgr A*, local mass/surface density, and the velocity dispersion in Baade's window, parameters of the potentials are 
brought up-to-date. The mass at large radii - in particular, that of the dark matter halo - is hereby constrained by requiring that the 
most extreme known halo, blue horizontal-branch star has to be bound to the Milky Way. 

Results. The Galactic mass models are tuned to yield a very good match to recent observations. The mass of the dark matter halo 
is - within the limitations of the applied models - estimated in a fully consistent way. As a first application, the trajectory of the 
hypervelocity star HE 0437-5439 is investigated again to check its suggested origin in the Large Magellanic Cloud (LMC). 
Conclusions. Despite their simplicity, the presented Milky Way mass models are very able to reproduce all observational constraints. 
Their analytical and simple form makes them ideally suited for fast and accurate orbit calculations. The LMC cannot be ruled out as 
HE 0437-5439's birthplace. 

Key words, dark matter - Galaxy: fundamental parameters - Galaxy: halo - Galaxy: kinematics and dynamics - Galaxy: structure - 
Stars: kinematics and dynamics 



1. Introduction 

Tracking large-scale orbits within the Milky Way is an impor- 
tant component in the study of the dynamical properties of stars, 
globular clusters, or satellite galaxies. Given the object's coor- 
dinates and velocity components, this is a straightforward task 
once the Galactic gravitational potential is available. Owing to 
a lack of observational constraints, a great deal of the Milky 
Way's mass distribution - in particular of its dark matter halo - 
remains unknown. Therefore, different mass models are able to 
reproduce the observations equally well. For the purpose of nu- 
merical orbit calculations, a mathematically simple and analyti- 
cally closed potential is preferred because it supports fast com- 
putations, which is especially favorable when using Monte Carlo 
methods for error estimation. 



The Galactic gravitational potential of Allen & Santillan 
( 1991| l perfectly fulfills these criteria and has thus been 
widely used to calculate trajectories for, e.g., globular clusters 



( |Odenkirchen et al p997l | Allen et al.||200^[Lane"et al.||2012), 
dwarf spheroidals dLepine e t al.||201l| , planetary nebulae (Wu 
|et al.|2011) , white dwarfs ([Pauli et al.||2003l |20"06"|), horizontal - 



branch stars (Altmann & de Boer 2000; Kaempf et al 
subluminous B stars (Altmann et al. 2004; Tillic h et al 



halo stars (|Schuster et a~ 
stars (Irrg ang et al.||2010[ 



2012 Pereira et al. 



Silva & Napiwotzki 2011 1, and hy 



2005 



2011 



2012 1, runaway 



pervelocity stars (Hirsch et al. 2005 Edelmann et al. 2005| l. 
However, the model parameters can be recalibrated using new 



and improved observational constraints. In this context, a sim- 
plified expression for the halo component is introduced as well. 

As for similar studies, standard constraints encompassing 
terminal velocities for the inner rotation curve, circular veloci- 
ties of maser sources for the outer rotation curve, in-plane proper 
motion of Sgr A*, local mass/surface density, and the velocity 
dispersion in Baade's window are fitted to determine the respec- 
tive parameters. To constrain the mass at large radii, in particu- 
lar that of the dark matter halo, the approach of Przybilla et al. 
( |20 1 0) > is followed here by requiring the dynamically peculiar 
halo star SDSSJ153935.67+023909.8 (J1539+0239 for short) to 
be gravitationally bound to the Milky Way (Sect.[2]l. 

Despite its simplicity, the resulting revised [Allen & Santillan 
model (hereafter denoted Model I) is capable of meeting all of 
its imposed conditions. In addition to the fast computation of 
realistic trajectories, calibration of the potential using the star 
J 1539+0239 allows the halo mass to be estimated in a fully con- 
sistent way (Sect.[3]l. To investigate the influence of the applied 
halo component on the results, the fitting is redone with the orig- 
inal halo mass distribution being replaced by a truncated, flat 
rotation curve model according to |Wilkinson & Evans| {1999 1 



(hereafter denoted Model II, Sect.|4|l and a |Navarro et al.| ( |1997 1 
dark matter halo (hereafter denoted Model III, Sect. [5). 

The spatial origin of the hypervelocity star HE 0437-5439 
(Edelmann et al. 2005) is investigated as a first application ex- 
ample (Sect.|6|. Finally, conclusions are drawn (Sect. |7]i. 
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2. Observations and fitting 

In this section, all the details of the fitting process are outlined, 
i.e., what assumptions are made, what observations are taken 
into consideration, and how the fitting actually is done. 



2.1. Introductory remarks 

All gravitational potentials 0(r, z) considered here are the sum 
of a central spherical bulge component <&b(K), an axisymmetric 
disk i>d(r,z), and a massive spherical dark matter halo Oh(^): 



0>(r,z) = fcbCRfoz)) + *d(r,z) + . 



(1) 



Here, (r,<p,z) are cylindrical coordinates and R = R(r,z) = 
Vr 2 + z 2 is the spherical radius. 

The corresponding density can be derived by virtue of 
Poisson's equation: 



4jrGp(r,z) = V 2 <D(r,z). 



(2) 



For the sake of simplicity, all potentials are independent of 
time. Consequently, rotating, nonaxisymmetric features, such as 
the Galactic bar or the spiral arms, are not implemented here, be- 
cause they require that more complex models (see, e.g., Pichardo 
|et al.|[2003j |2004| l are considered. To account for the effects 
of bar and spiral arms, additional systematic uncertainties are 
added to the observational data where necessary, see for instance 
Sect. l232l 

Sticking to the convention of Allen & Santillan ( 1991| l, 
throughout this paper, the gravitational potentials are expressed 
in units of 100km 2 s -2 , lengths in kpc, and masses in galactic 
mass units M ga i = 100 x 1000 2 kpc [m]/G [SI] M© * 2.325 x 
10 7 Mq, yielding a gravitational constant G of unity. 

2.2. Observational constraints 



Similar to |Dehnen & Binney ( 1998 1, observations are divided 
into eight different groups that are introduced and discussed in 
the following. 

2.2.1. Solar kinematics 

Analyzing kinematic observations made in the celestial refer- 
ence frame in a Galactocentric picture requires a coordinate 
transformation from one system to the other. Since this trans- 
formation obviously depends on the Sun's position and veloc- 
ity, the interpretation of any data is affected by the values cho- 
sen for these parameters. To set up the coordinate transforma- 
tion, we take J2000.0 coordinates for the Galactic center (GC) 
a GC = 17 h 45 m 37 s .224, 6 GC = -28°56' 10723 and north Galactic 
pole (NGP) a NG P = 12 h 51 m 26 s 282, 6 NGP = +27°07 , 42 , . , 1, as 
mentioned in the appendix of Reid & Brunthaler ( 2004| l. By 
monitoring stellar orbits around the central supermassive black 
hole Sgr A*, Gillessen et al. (2009i restricted the distance r© 
from the Sun to the GC to the range 



r© = 8.33 + 0.35 kpc, 



(3) 



which is the first observational constraint. 

The Sun's velocity is a superposition of the in-plane cir- 
cular motion vo of its local standard of rest (LSR) around the 
GC and its peculiar motion vq = (U, V, W)q relative to the 
LSR. Hereby, U is the component towards the GC, V in direc- 
tion of Galactic rotation, and W perpendicular to the Galactic 



plane. Motivated by the McMill an & Binney] ( |2010| l discus- 
sion of a systematic offset in the V-component of maser mo- 
tions versus a nonstandard value of V©, the (U, V, W)q values 
determined by Sch onrich et al.| ( |2010| l are in this case prefer- 
able to others since they allow the use of maser sources to con- 
strain the Galactic rotation curve (see Sect. |2.2.2| i without the 
necessity of modeling an additional net peculiar motion. Thus, 

v = (H-l!o75, 12 - 24 -o:47> 7 - 25 -o36) ± 2, 0.5) km s" 1 is used 
here. 



Following the argument in McMillan (201 1 1, the proper mo- 
tion of Sgr A* along Galactic longitude //sgrA* = -6.379 ± 
0.026 mas yr~' ( Reid & Brunthaler||2004 1 is utilized as an indi- 



cator for vo instead of Oort's const ants A and B. With a spatia l 
motion of Sgr A* of ± 1 kms" 1 ( jMcMillan & Binney||2010| , 
this proper motion can be assumed to be solely caused by the 
solar motion around the GC yielding as second observational 
constraint: 



MsgrA* = — - (Vq + vo) = — - (Vq + v c (r©)) 
rQ r Q 



(4) 



The circular velocity v c (r) is hereby linked to the potential via 



,d<D(r,0) 



dr 



(5) 



2.2.2. Rotation curve 

The Galactic rotation curve deduced from terminal velocities in 
the inner and intermediate regions and from maser sources in the 
outer region yield the third, fourth, and fifth groups of observa- 
tional constraints considered here. 



Terminal velocities: the terminal velocity v term (l) is the mea- 
sured peak velocity of the interstellar gas along the line-of-sight 
at Galactic coordinates b — and -90° < / < 90°. Assuming a 
circularly rotating interstellar medium, v term (/) can be computed 
from the circular velocity v c at r — r© sin(Z) corrected for the 
observer's projected motion: 



v t ei-m(0 = v c (r© sin I) - (v c (r©) - V©) sin / - Uq cos / . 



(6) 



For this study, unprocessed terminal velocities are extracted 
from surveys in H i by Burton & Gordon ( 1978[) and in CO by 
|Clemens| ( fl985| l. Similar to |Dehnen & Binney| ( |1998] l, a constant 
uncertainty of 7 km s" 1 is added in quadrature to the stated uncer- 
tainties of the IClemensI data to account for non circular motions 
and the inability of the model potential to reproduce, e.g., spiral 
arm features. Complementary to the CO measurements, the Hi 
data probe the central region with | sin l\ < 0.3, which is addi- 
tionally affected by distortions due to the bar. Therefore, a rather 
generous uncertainty of 15 km s _1 is quadratically added to their 
measurement errors. 



Maser observations: maser sources are associated with high- 
mass star-forming regions and are consequently good tracers 
of the kinematic properties of their surrounding interstellar gas. 
Using very long baseline interferometry techniques, their posi- 
tions, parallaxes, and proper motions can be measured to very 
high accuracy even for distant sources giving the opportunity to 
sample a relatively wide range of Galactocentric radii. Together 
with heliocentric radial velocities from Doppler shifts, these 
data hence allow precise determination of the masers' three- 
dimensional motions. In contrast to Reid et al. ( 2009 ), McMillan 
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|& Binney| ( |2010| l show that it is fairly justified to adopt their 
corresponding circular velocity components as probes for the 
Galactic rotation curve when taking the Schonrich et al. val- 
ues for the peculiar motion of the Sun. Following this argu- 
ment, parallaxes, proper motions, and heliocentric radial veloci- 
ties for 30 maser sources have been compiled from the literature 
USandstrom et al .pMf\ |Hirota et al.|2008a|b| [Reid et al.|2009| 
Rygl et al.||2010[ |Sato et al.||2010| |Asaki et al.||2010[ |Niinuma] 



et al.|20iTJ Honma et al.|201 To account for the virial motion 



of the masers with respect to their high-m ass star-forming r e 
gions, a constant uncertainty of 7 km s (see Reid et al. 2009 [ ) is 
added in quadrature to the measurement errors when computing 
circular velocities. 



2.2.3. Mass and surface densities 

The local dynamical mass density pq and the surface den- 
sity Si i are together the sixth group of observational constraints. 
Their derivation is straightforward from the model: 



PQ = Pb(r©) + Pd(r©, 0) + p h (r©) 

1.1 kpc 

Ei.i = J (Pb(rQ,z) + Pd(rQ,z) + Ph(rQ,z))dz. 



(7) 



(8) 



-1.1 kpc 



Holmberg & Flynn[ obtained p = 0.102 ± 0.010 M Q pc 3 from 
HrppARCos data on a volum e-complete sample of A and F stars 
( |Holmberg & Flynn|2000"| and £u = 74 + 6M pc~ 2 from ob- 
servations of K giant stars ( Holmberg & Flynnp 004 ). 



2.2.4. Velocity dispersion in Baade's window 

The velocity dispersion of the bulge in Baade's window, <x BW , 
is implemented as the seventh constraint probing the inner most 
region of the Milky Way. Following the reasoning of Dehnen & 
|Binney| ( fl998| l, it is given as 



cr BW = cr b (O.O175r , -0.068r o ) =117+ 15kms~ 
where cr b is estimated from the model potential as 

oo 

°b( r .z)= — : — r Pb(r,z) — — — az. 
Pb(r,z) J dz 



(9) 



(10) 



2.2.5. The extreme halo star J1 539+0239 

The most stringent constraints on the halo mass result from the 
kinematically most extreme objects. Usually, the very distant 
satellite galaxy Leo I is considered to be the most important 
tracer of the halo mass. Its proper motion has recently been 
measured with the Hubble Space Telescope (HST) for the first 
time (So hn et al.|2012 1 and allowed to determine a limit on the 
halo mass (Boylan-Kolc hin et al.|2012) . However, close-by halo 
stars may provide even more compelling results. The kinemat- 
ically most extreme halo star is the blue horizontal-branch star 
J1539+0239 for which |Przybilla et al.| < |20T0| > found a Galactic 
rest-frame velocity of about 69A + j^ k m s~\ making it the fastest 
halo object known so far. Unlike other hypervelocity stars, how- 
ever, this object is currently approaching the Galactic disk with 
its pericentric passage of the GC still to occur. As discussed in 
the next paragraph, this indicates that J1539+0239 is gravitation- 
ally bound to the Milky Way, yielding a so far unexploited and 



- due to the star's immense space motion - significant constraint 
on the Galaxy's total mass. We show in Sects. [3~2 4.2 and 5.2 
that the limits on the halo mass derived from the kinematics of 
this star alone are even more stringent than those from Leo I, 
so we add J1539+0239 to our set of constraints for the poten- 
tials and use Leo I as a consistency check. Our approach, which 
is based on the six-dimensional phase space information of the 
most extreme objects, is complementary to statistical studies of 
several thousand stars (see, e.g., Xue et aL]|2008| |Gnedin et al. 
2010) that only make use of line-of-sight velocities. 

To justify the assumption that J 1539+0239 is bound to the 
Galaxy, consider the opposite, i.e., a trajectory that passes by 
the Milky Way only once, never to come back again. Such an 
orbit could be explained either by a non-Galactic origin, e.g., 
as a hypervelocity star ejected from another galaxy, or by an ex- 
treme dynamical event, e.g., a very close multibody encounter in 
a globular cluster or satellite galaxy, ejecting a star from its halo 
environment in a direction pointing not too far away from the 
GC. Although both scenarios cannot be ruled out by stellar age 
or space motion arguments, they are unlikely since they require 
very special circumstances, e.g., an extraordinarily high ejection 
velocity, a specific ejection direction, or the ejection event occur- 
ring at a certain time. Most of these restrictions are weaker or do 
not apply at all in the case of a bound orbit. For instance, the pe- 
culiar motion of J1539+0239 could still be the result of an ejec- 
tion event. But because the accessible phase space is finite now, 
there is a non-zero probability of finding a trajectory starting 
from any initial condition (with matching conserved quantities) 
sooner or later in a phase-space state very close to J1539+0239, 
allowing for countless more ejection locations, directions and 
times. In particular, the GC no longer can be excluded as the 
star's spatial origin, which opens up the possibility of a dynam- 
ical interaction with the central supermassive black hole (Hills 
1988). This mechanism is known to be most powerful in terms 
of acceleration and is the best candidate for explaining extreme 

( |2009> 



velocity objects. Alternatively, as shown by Abadi et al 



the latter may result from the tidal disruption of a dwarf galaxy, a 
non ejection scenario that is also able to produce bound (but not 
unbound), high-velocity stars that approach the Galactic disk. 

The assumption that J1539+0239 is gravitationally bound to 
the Milky Way is the eighth and last observational constraint. 
With E t denoting the star's total energy, i.e., the sum of kinetic 
and potential energy, per unit mass, the standard approach to in- 
corporating this is to require E„ < 0. But the Milky Way is not 
an isolated system, and beyond some distance, the influence of 
other objects such as Andromeda is no longer negligible. Since 
those effects are not included in the presented models, a more 
rigorous approach is chosen here by defining a bound state as 



E, < (D(0,200kpc), 



(11) 



theoretically limiting orbits to R ~ 200 kpc. However, the actual 
implementation in the least-squares fit is 



£*- O(0, 200 kpc) 



X* = 



AE, 







if E, - 0>(0,200 kpc) > 0, (12a) 
otherwise. (12b) 



Hereby, AE, is the error in E t resulting from uncertainties 
in determining J1539+0239's distance, radial velocity, and 
proper motions. Since AE„ is large (for instance E*/AE* ~ 
0.82 in the original Allen & Santillan potential), violations of 
Eq. ( fTT) i are generously admitted. For the same reason, the 
particular choice of 200 kpc is not decisive as for instance 
-O(0,200kpc)/A£„ ~ 0.16 while -O(0, 150kpc)/A£„ ~ 0.21 
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and -D(0 ,250kpc)/A£ , < , ~ 0.13 (again in the original |Alien~&l 



Santillan potential). 



2.3. Fitting process 

The model parameters (see Sects. [3jj5]) are simultaneously de- 
termined via a ^-minimization procedure using the Interactive 
Spectral Interpretation System (ISIS, Houck & Denicola 2000). 
Observations of the rotation curve, as well as of J1539+0239, are 
always taken in raw form from the literature, i.e., as measured in 
celestial coordinates, in order to exclude inconsistencies during 
the transformation to the Galactic reference frame originating in 
different values for the LSR and solar motion. Moreover, any 
given uncertainty found in the literature is considered here, typi- 
cally via Gaussian statistics. For instance, to compare the maser 
measurements to model predictions, their positions and veloci- 
ties are converted to Galactic coordinates, allowing computation 
of rotational velocities v c and Galactrocentric distances r. The 
corresponding uncertainties in the observational data are then 
accounted for using Gaussian error propagation. Additionally, 
the (admittedly small) uncertainties stemming from vq are prop- 
agated in the same way. As the resulting data points in the 
(r, v c )-plane have errors in both directions (Ar, Av c ), it is not 
obvious at which radius to compare the model to observation, 
i.e., how to compute the x 2 without neglecting Ar. The solu- 
tion chosen here is to convert Ar to an error in v c by estimat- 
ing the effects of it on the model rotation curve according to 
(Av c ) mol jei = (dv c /dr) mo deiAr and adding it in quadrature to Av c . 

Very much attention is paid to take all sources of error into 
consideration in order to assign the correct weighting to each ob- 
servational constraint in the total x 2 ■ Unfortunately, systematic 
effects such as the missing modeling of spiral arms are difficult 
to quantify and are not sufficiently dealt with by just adding a 
constant error. Therefore, an unweighted fitting routine, e.g., as- 
suming that a single maser point weighs the same as the proper 
motion of Sgr A* (see Eq. Q), does not yield a satisfactory re- 
sult because the fit would be dominated by the large number of 
data points used while fitting the rotation curve. To account for 
this, each contribution to the total x 1 of the above-mentioned 
eight groups of observational constraints is divided by the num- 
ber of data points N within the group: 



8 ( 1 Ni 



(13) 



The model parameters given in Tables [TJ|3] are those that mini- 
mize this weighted^- 2 . 

3. Model I 



Model I is a revision of the Allen & Santillan ( 1 99 1 ] > potential 



3. 1 . Components and characteristics 
3.1.1. Bulge and disk 

The potential forms <t>b(R) and 3>d(r, z) of bulge and disk are 
those proposed by Miyamoto & Nagai ( 1975 1: 



fb(fi) = - 



V« 2 + h 2 



fd(r,z) = - 



^r 2 + (a d + Vz 2 + ^d 2 ) 



(15) 



The parameters M\,/M d determine the contribution of their com- 
ponents to the total potential, and b\,/a d /b d are scale lengths. 
The related densities pb and pd are (see Eq. Q) 



p b (R) = 



Pd(r,z) 



3b h 2 M b 



An(R 2 + b h 2 ) 5 ' 2 



2 M d a Ar 2 + (fl d + 3 Vr 2 + z 2 ) (a d + Vr 2 + z 2 f 



An 



(z 2 + b d f 2 (r 2 + (a d+ ^^d 2 )) 



2\5/2 



(16) 



(17) 



Integrating these densities over the entire volume to obtain the 
total masses and m d of bulge and disk gives the expected 
identities = Mb and m d = M d . 

3.1 .2. The dark matter halo 

Pursuing mathematical simplicity, the halo potential <t>h(R) used 
here differs slightly from the original one by |Allen "& Santillan 
We assume the halo mass «h inside a sphere of radius R is given 
by (see |Allen & Martos|1986} 



m h (< R) = 



i HIT 1 



it) 



if/?<A, (18a) 



= constant, otherwise. (18b) 



Here, is again a weighting factor and a scale length. 
Clearly, this specific form is chosen to have the asymptotic be- 
havior mh(< R) cc R for large R motivated by the observed flat 
rotation curve. To avoid an unphysical, infinite halo mass, a cut- 
off parameter A is incorporated as well. In principle, the expo- 
nent y is a free parameter, too. However, the shortage of halo 
constraints means the ability of the model to reproduce the ob- 
servations is not attenuated by setting y = 2, reducing the ex- 
pression's complexity (and avoiding singularities at the origin in 
the equations of motion for y < 2). 

The expression nih(< R) is the volume integral of the halo 
density ph over the sphere of radius R. Assuming a spherically 
symmetric halo density ph = Ph(^) and using Poisson's equation 
(Eq. Q), a relation between »ih(< R) and its respective potential 
can be derived: 



R R 

m h (<R) = An j R 2 p h (R)dR = j R 2 V 2 ii> h <RhlR 

o o 



= f R 2 ^A\R 2 A®h(R)\dR 
J R 2 dR\ dR ) 



dR h ' 



(19) 



Mzn Inserting Eq. ( 18 i and directly integrating the result accounting 
for the boundary condition Oh(°o) = yields 
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«h (kpc) 

Fig. 3. Visualization of the correlation between Mj, and ah in 
Model I: The single parameter 90%-confidence region defined 
by Ax 2 < 2.71 and corresponding to a 74%-joint-confidence re- 
gion is black-rimmed and shaded in dark gray, while the 3cr re- 
gion (Ax 2 < 9, 99%-joint-confidence region) is the total black- 
rimmed, gray-shaded area. The cross marks the location of the 
best fit. The four curves define the loci of constant halo mass 
inside a sphere of radius 14. 1 kpc (red dashed), 50 kpc (blue dot- 
ted), 100 kpc (blue dashed-dotted), and 200 kpc (blue dashed- 
dotted-dotted). 



A4 



1 



(r-i) 



In 



f'+(£) 



y-l \ 



1 + 



1 + 



(if 



^ + ( a h ) 



y-l ' 



ifR < A{20a) 



otherwis(20b) 



For A = 100 kpc and y — 2.02, this is equivalent to the expres- 
sion given in |Allen & Santillari ( |1991| l obtained via integration 
by parts. The advantage of the above representation is that the 
dependent variable R appears only once and linearly in the argu- 
ment of the logarithm (since y = 2). 
The corresponding density ph is 



Ph(R) 



(sr'ftsr+r) 



R 2 



ifR < A, 



otherwise. 



(21a) 



(21b) 



3.2. Results 



The properties of the best-fitting Model I are visualized in Figs.[T] 
and|2]and summarized in Table[T] Given the model's simplicity, 
the overall agreement with the imposed constraints is very good. 
In particular, rQ, jUsgrAs PQ, £1.1, and <x BW are almost perfectly 
reproduced. The rotation curve, on the other hand, reveals sys- 
tematic discrepancies between model and observation originat- 
ing mainly in the inadequate treatment of the Galactic bar and 



the spiral arms. Moreover, most of the maser sources still seem 
to lag behind the Galactic rotation curve. Despite these short- 
comings, the unweighted x 2 P er degree of freedom (d.o.f.) is 
about 1.5, thereby confirming the good match. 

To see how well the individual parameters are confined by 
the observational constraints, 90%-confidence limits were com- 
puted from the weighted x% statistics after multiplying it with 
a factor yielding ^/d.o.f. = 1 at the best fit, i.e., at the min- 
imum^. The resulting intervals are given in Table [I] too. As 
illustrated in Fig. 3] there is a strong correlation between Mh and 
Oh that explains the large uncertainties of these two parameters. 
This degeneracy is a consequence of the missing constraints on 
the halo mass distribution. The available observations are only 
able to constrain the halo mass inside a sphere of radius 14.1 kpc 
(see Fig. [3]), which, however, can be reproduced by numerous 
different combinations of M^-a^ pairs. A parametrization of the 
halo using only one independent variable would thus be suffi- 
cient to obtain a good fit to the data. From Fig. [3] it is also clear 
that the degeneracy in the two halo parameters does not follow 
curves of constant halo mass inside 50, 100, or 200 kpc, which 
explains the relatively large uncertainties of the total masses in- 
side these radii as listed in Table Q] 

Przybilla et al. ( 20 1 0] > were the first to exploit the star 



J 1539+0239 as a probe for the mass of the Galaxy and its 
main contributor, the dark matter halo, and they derived mh > 
1.7^j j x 10 i2 Mq based on the same model potential as dis- 
cussed here. Their result, however, was derived solely from the 
condition that J1539+0239 is bound to the Milky Way while ne- 
glecting any other observational constraint that keeps the mass 
low. In contrast, the mass estimates presented in Table[T]are fully 
consistent with all the observational data mentioned before. 

The best-fit parameters yield a total mass (with 90%- 
confidence limits computed analogously to the model parame- 
ters' uncertainties) of M tota] = 1. 9* n g X 1O 12 M0 « m h . This 
result is a modest revision of the Przybilla et al.| ( |2~010| l value, 
and it confirms their conclusion that a total Galactic mass of 



about 1 x 10 i2 Mq such as found by Xue et al. ( 2008 1 possi- 
bly underestimates the true value. From Eq. (|18|l and the fact 
that A = 200 kpc at the best fit, it is obvious that M tota i is af- 
fected by the restriction A < 200 kpc. Although the latter is 
justified by cosmological simulations on galaxy formation, the 
specific choice for the upper boundary can of course be subject 
to discussion. But since even the mass corresponding to A = 
450 kpc lies within the 90%-confidence interval of M tota i, this 
concern does not seriously change the result. To remove the A 
dependence, consider the quantity M«<50k P c, i-e-, the total mass 
within 50 kpc. From the kinematic analysis of a combined sam- 
ple of field horizontal-branch stars, globular clusters, and satel- 
lite galaxies, Sakamoto et al. ( 2003} have obtained Mfl<5oicpc = 
5.4^ x kWm© agreeing with M fi<5 okpc = 5.1^.3 x 10 « Mq 
as found in this study. 

It is worth noting that J 1539+0239 is not at all gravitation- 
ally bound to the Milky Way in the best-fitting model. This is 
possible because of the large uncertainty in its position and ve- 
locity determination resulting in a small contribution (x* = 0.66) 
to the total x\- More precise kinematic information on the star 
would set tighter constraints on the Galactic potential, probably 
reducing the uncertainties in the mass estimates. 

We also model the trajectory of the very distant satellite 
galaxy Leo I based on the full six-dimensional phase space in- 
formation (position: a = 10 h 8 ra 28^68, 6 = +12° 18' 19.7"; dis- 
tance: 256.7 ± 13.3 kpc; radial velocity: 282.9 ± 0.5kms _1 ; 
proper motion: fi a cos 6 = — 0. 1 140 + 0.0295 mas yr 1 , fig = 



-0.1256 ± 0.0293 mas yr 1 , Sohn et al. 2012i. Leo I is found 
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Table 1. Properties of the best-fitting Model I. 



Parameter 



Value 



Constraint 



Value 



Observation 



Model 



Derived Quantity 



Value 



r Q (kpc) 


8.40 ± 0.08 


M b (M gal ) 


409 ± 63 


M d (M gal ) 


2856 + ^Q2 


M h (M gal ) 


1018^™ 3fl 


b b (kpc) 


0.23 ± 0.03 


a d (kpc) 


A 99+0.53 
^•^ -0.99 


b d (kpc) 


n 909+0.020 

"'"^-0.025 


a h (kpc) 


9 c^9 +25.963n 
iJui -1.4U 


A (kpc) 


200+°/ 


y 


2 (fixed) 



r© (kpc) 
/"SgrA* (masyr -1 ) 
Terminal velocities from H i 
Terminal velocities from CO 
Circular velocities from masers 
p (Mo pc- 3 ),E M (M© pc- 2 ) 
cr BW (kms"') 
X* 



8.33 + 
-6.379 : 



0.35 
: 0.026 



8.40 
-6.384 



See Fig. 
See Fig. 
See Fig. 
0.102 ±0.010, 74: 
117 ± 15 
<0 



See Fig. 
See Fig. 
See Fig. 
0.102, 74 

120 

0.66 



v (kms ') 
mb (10 9 M o ) 
md (10'°M o ) 
mh (10 12 M o ) 
M R<5 okpc(10 12 M ) 
M R<10 okpc(10 12 M ) 
M s<2001cpc (10 12 Mo) 
v csc (kms -1 ) 
A (km s -1 kpc -1 ) 
fi(kms -1 kpc -1 ) 



242.0 
9.5 ± 1.5 

6 6 +0 - 9 



1. 



3+2-4 

51 +0 - 33 

"• J1 -0.04 
A Q7+0.96 
"• 7/ -0.09 

1 9 +2 - 4 
616.4 
15.06 
-13.74 



Notes. The quoted uncertainties for the model parameters and masses are 90%-confidence limits: After normalizing the weighted^ 2 via multi- 
plication with a factor yielding,^ 2 /d.o.f. = 1 at the minimum^ 2 , 90%-confidence intervals are calculated from this normalized^- 2 statistics and 
the condition Ay 2 = 2.71. (<l) The large uncertainties are due to a strong correlation between M h and a h , see Fig. [5] (W Motivated by cosmological 
studies, the restriction A < 200 kpc is imposed. 
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Fig. 1. Comparison of the best-fitting model rotation curve (solid black line) with terminal velocities from surveys in H i (red circles) 
and in CO (green diamonds), as well as maser observations (blue squares) for Model I. 



to be formally bound to the Galaxy, showing that J1539+0239 is 
indeed the dominant constraint on the halo mass. 



Table [T also contains some derived quantities that are not di- 
rectly fitted, such as v esCj Q or vq/tq. The escape velocity at the 
Sun's position is about 616kms _1 , thereby lying very close to 
the 90%-co nfidenc e limit 498 km s 1 < v esc < 608kms _1 of 
(2007). The ratio of vo over rQ is 28.8 km s -1 kpc" 



Smith et al. 



and th us "somewhat less than the values found by |Reid et al 
( |2009| l, 30.3 ± 0.9 km s 1 kpc -1 , or |McMillan & BinneyfpOTUl r 



29.9 - 3 1 .6 km s kpc , from maser sources alone. 



As pointed out before, the halo mass distribution is not well 
constrained by the observations, so different potential shapes of 
the halo can yield almost equal matches to the data but differ 
significantly in other properties, as for instance the total mass. To 
investigate this behavior, two other representations of the dark 
matter halo are considered in the following sections. 



4. Model II 

In Modell II, the halo is replaced by the truncated, flat rotation 
curve model used in Wilkinson & Evans (1999) and Sakamoto 
[itaLlpOOl) . 

4. 1 . Components and characteristics 

While bulge and disk components have the same shape as in 
Model I, the halo potential reads as (Wilki nson & Evans[i 999) 



O h 0R) 



ln 



f V# 



R 



The resulting density 
M h 



p h (R) = 



a h - 



4^ & (R2 + flh 2) 



2^/2 



(22) 



(23) 



is cusped like R 2 for R <K ah and falls off like R~ 5 for R » 
ah- In this way, the corresponding rotation curve is flat in the 
inner regions, and the total halo mass nih = Mh is finite without 
invoking a cutoff parameter. 
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Fig. 2. Gravitational potential <b(r, z), formal escape velocity v esc (r, z) = V-2<t>, and total mass density p(r, z) for the best-fit param- 
eters of Model I (upper panel), Model II (middle panel), and Model III (lower panel). Contours are indicated by dashed lines. 



4.2. Results 

The properties of the best-fitting Model II are visualized in 
Figs. [2] and [4] and summarized in Table [2] While the rotation 
curve reveals the same systematic shortcomings as found in 
Model I, the match of the remaining constraints is of comparable 
quality. The unweighted^ 2 per d.o.f. is 1.7. 

Again, there is a tight correlation between the halo param- 
eters Mh and Oh (see Fig. [5]). To avoid unphysically large ha- 
los, the restriction a^ < 200 kpc has to be imposed based on 
exactly the same reasoning as the condition A < 200 kpc in 
Model I. Owing to the correlation, this step also sets an up- 
per limit to the parameter Mh, hence to the total halo mass 
nib. Unfortunately, the effects on the latter and thus on the to- 
tal Galactic mass are more severe in this case since its 90%- 



confidence interval, M tota i = 1.7^2-? x 10 12 Mq, covers a^ val- 
ues only up to 218 kpc. The mass inside 50kpc, on the other 
hand, is almost independent of the chosen upper boundary for 
«h- According to Sakamoto et al. (2003), this quantity is much 
more robust than, e.g., the total mass of the Milky Way. The 
reason for this is the specific form of the applied potential. It 
is such that loci of equal masses inside spheres of radii rang- 
ing from ~0 to ~50kpc lie in a very narrow band in the Mj,- 
«h plane (see Fig. [5}. Although the correlation follows contours 
of equal mass within the central ~20kpc, it therefore nearly 
follows contours of equal mass within the central 50 kpc, re- 
sulting in small uncertainties for the latter. The value derived 
here, M R<50kpc = 4.6 + 0.3 x 1O U M0, is slightly lower than in 



Sakamoto et al. 



(2003 1, M R<50kpc =5.4^x10" M, 



i . 
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Table 2. Properties of the best-fitting Model II. 



Parameter 



Value 



Constraint 



Value 



Observation 



Model 



Derived Quantity 



Value 



r Q (kpc) 


8.35 ± 0.08 


M b (M gal ) 


175 ± 28 


M d (M gal ) 


2829 ± 192 


M h (M gal ) 


69 725+5790 a 
U7 / —20 93 1 


b b (kpc) 


0.184 ±0.040 


a A (kpc) 


4 85 +0 41 

t'° -0.33 


bi (kpc) 


0.305 ± 0.020 


a h (kpc) 


200+° fl - h 



-') 



r© (kpc) 
yU SgrA . (mas yr~ 
Terminal velocities from H i 
Terminal velocities from CO 
Circular velocities from masers 
p (Mo pc- 3 ),E M (M Q pc- 2 ) 
cr BW (km s~') 
X. 



8.33 ±0.35 
-6.379 ± 0.026 



8.35 
-6.383 



See Fig. 
See Fig. 
See Fig. 
0.102 ±0.010, 74: 
1 17 ± 15 
<0 



See Fig 
See Fig 
See Fig 
0.102, 75 

116 

0.80 



v (kms ') 
m„(10 9 M o ) 
md (10'°M o ) 
mh (10 12 M o ) 
M R<5 okpc(lO 12 M ) 
M R<10 okpc(10 12 M ) 
M« <2 ookpc(10 12 Mo) 
v esc (kms^ 1 ) 
A (km s -1 kpc -1 ) 
S(kms-'kpc _1 ) 



240.4 
4.1 ±0.7 
6.6 ±0.5 

1 6 +0 2 

1 - u -0.5 

0.46 ± 0.03 

n 70+O.O6 
"• /7 -0.08 

1 2 +0 1 
575.9 
15.11 
-13.68 



Notes. The quoted uncertainties for the model parameters and masses are 90%-confidence limits (see notes on Table[T]for details). {a) The large 
uncertainties are due to a strong correlation between M b and a b , see Fig. [5] {b) Motivated by cosmological studies, the restriction aj, < 200 kpc is 
imposed. 




10 



100 



r (kpc) 



Fig. 4. Comparison of the best-fitting model rotation curve (solid black line) with terminal velocities from surveys in H 1 (red circles) 
and in CO (green diamonds), as well as maser observations (blue squares) for Model II. 



In general, the gravitational potential of Model II is shal- 
lower than in Model I (see Fig. [2]) implying systematically lower 
masses, a lower local escape velocity of about 576 km s -1 in 



agreement with Smith et al. ( 2007| ), and an unbound orbit for 
J1539+0239 (x* = 0.80). Nevertheless, Leo I is still on a bound 
orbit. 



5. Model III 

The third halo potential considered here is based on the universal 
density profile of dark matter halos suggested by Navarro et al. 
1997| i from cosmological simulations. 



5. 1 . Components and characteristics 

Using 



Ph(R) = 



M h 



1 



4?r (a h +R) 2 R 



(24) 



(Nav arro et al.|[l997 i together with Eq. d2J, one obtains a halo 
potential of 



O h 0R) 



M h , /. R 

In 1 + — 

R \ a h 



(25) 



The weighting factor is thereby equivalent to the mass inside 
a sphere of radius ~5.3 times the scale length ah- The combi- 
nation of this halo potential with a Miyamoto & Nagai (1975 ) 



disk and bulge component (see Section 3.1.1 1 is hereafter de 



noted Model III. Owing to the use of Eq. (|24|l for R — > 00, the 
total halo mass is formally logarithmically divergent. 



5.2. Results 

The properties of the best-fitting Model III are visualized in 
Figs. [2] and [6] and summarized in Table [3] The quality of the 
fit is similar to Models I and II and yields an unweighted x 1 per 
d.o.f.ofl.7. 

In contrast to the two previous models, the rotation curve is 
rising outside the solar circle, reaching its maximum at about 
82 kpc, and falling off beyond this point. This is a consequence 
of the search for a compromise between the opposing constraints 
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Table 3. Properties of the best-fitting Model III. 



Parameter 



Value 



Constraint 



Value 



Observation 



Model 



Derived Quantity 



Value 



r Q (kpc) 
M b (M gal ) 
M d (M gal ) 
M h (M gal ) 
b b (kpc) 
a& (kpc) 
b A (kpc) 
a h (kpc) 



8.33 ± 0.09 



439: 

3096: 



28 
197 



142 200 +137900a 
0.236 ±0.021 

-t 9AO+0.144 
''■ iui -O.I21 

0.289 ± 0.022 

,J,U -16.78 



r e (kpc) 
fi SgrA , (masyr" 1 ) 
Terminal velocities from H i 
Terminal velocities from CO 
Circular velocities from masers 
p© (M© pc" 3 ), Si,! (M© pc" 2 ) 
o- BW (kms" 1 ) 
X* 



8.33 ± 
-6.379 : 



0.35 
: 0.026 



See Fig. I 
See Fig. i 
See Fig. I 
0.102 ±0.010, 74: 
117 ± 15 
< 



8.33 
-6.380 



See Fig 
See Fig 
See Fig 
0.102, 75 
123 
0.20 



239.7 
10.2 ±0.7 
7.2 ±0.5 



v (kms ') 

m b (lO 9 M0) 

m d (10 10 M o ) 

m h (10 12 M o ) 
M s<50kpc (10 12 M o ) 
M R<100kpc (10 12 M o ) 1.67 ±0.46 
M^ookpc d0 12 M o ) 

v e S c.O (km s "') 
A (kms -1 kpc" 1 ) 
B(kms -1 kpc" 1 ) 



0.81!°;! 3 



3 +1 2 
811.5 
14.70 
-14.08 



Notes. The quoted uncertainties for the model parameters and masses are 90%-confidence limits (see notes on Table[T]for details). {a) The large 
uncertainties are due to a strong correlation between M b and Oh, see Fig. [7] {b) Formal divergence, see Sect. 5.1 
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Fig. 6. Comparison of the best-fitting model rotation curve (solid black line) with terminal velocities from surveys in H i (red circles) 
and in CO (green diamonds), as well as maser observations (blue squares) for Model III. 



of making J1539+0239 bound, i.e., giving the Milky Way a very 
high total mass, and, at the same time, limiting the mass in the 
central ~ 19 kpc (see Fig. [7]) to match the remaining observa- 
tions. Unlike Models III, I, and II were designed to give flat 
rotation curves, avoiding any analogical behavior. As a result, 
the masses of Model III systematically exceed those of the oth- 
ers, e.g., by more than a factor of 1.5 in terms of M^gokpc- 
Moreover, because of using Eq. (24i for R — > oo, the total 



mass Mtotai is logarithmically infinite. Following Navarro et al. 
( |1997 1, an alternative measure of the halo mass is given by 
M200, which is the halo mass inside a sphere of radius r2oo de- 
fined to have a mean interior density of 200 times the critical 
value for closure p cl -; t = 3H 2 /%7rG. Assuming a Hubble con- 
stant H of 73 km s _1 Mpc _1 , values of r2oo = 319+^ kpc and 
M200 = 4.0+} g x 10 12 M© are derived from Model III agree- 
ing with r 2 oo = 301 kpc and M 2 oo = 3.4 x 1O 12 M0 as used 
by Abadi etaLl ( |2009) . The respective local escape velo city of 
8 12 kms" 1 significantly exceeds the Smith et al. (2007i value 
but is comparable to the one in |Abadi et al.| ( |2009| l. Although 
this is high enough to keep Leo I formally bound, it does not 
hold J1539+0239 within 200 kpc fa = 0.20). 



Kenyon et al. ( 2008} considered a Galactic potential that 



only differs from Model III in the bulge component: the ex- 



pression R + b\, is used instead of V/? 2 + b h 2 . Putting special 
emphasis on the very central region of the Milky Way, the 
parameters of their model were tweaked to provide a reason- 
able match to observations from R = 5pc to R = 100 kpc. 
Compared to Model III, their resulting model is less massive 
(M200 = 1.6 x 1O 12 M0), which is also reflected by lower values 
for the local escape velocity (v esc = 635 km s" 1 ), the local mass 
and surface density (p© = 0.046 Mope" 3 , S u = 46M0pc~ 2 ), 
and the circular motion of the LSR (vq = 220 kms" 1 ). 



6. The hypervelocity star HE 0437-5439 revisited 

Given the smoothness of the presented potentials, a simple 
fourth-order Runge-Kutta method with adaptive stepsize con- 
trol is sufficient to reliably and efficiently solve the equations 
of motion (as given in Appendix [A} numerically. Various trial 
trajectories, including those of the LSR, were successfully cal- 
culated in order to test the self-written ISIS-implementation of 
the adopted integrator. The hypervelocity star HE 0437-5439 
( Edelmann et al.|2005 1 is analyzed as a first and very interesting 
application because the spatial origin of this star is still under 
debate. 
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Table 5. Results of the kinematic investigation of HE 0437-5439. 



ku n 




ivioaei i 






IVlOQcl 11 






A/TriH^I TTT 
IVlOQcl 111 




d p (kpc) 


7/ p (Myr) 


v p (kms 4 ) 


d p (kpc) 


7/ p (Myr) 


v p (km s _1 ) 


d p (kpc) 


Tp (Myr) 


v p (km s _I ) 


#la 


13.2 ±4.5 


-22+13 


640 ± 70 


13.2 ±4.5 


-22 ± 13 


640 ± 70 


13.1 ±4.4 


-22 ± 13 


640 ± 70 


#lb 


15.6 ±10.4 


-98 ± 15 


690 ± 50 


15.7 ± 10.3 


-98 ± 16 


680 ± 50 


14.8 ± 10.1 


-95 ± 14 


730 ± 40 


#2 


11.0 ±4.5 


-23 ±12 


680 ± 80 


10.9 ±4.5 


-23 ± 12 


680 ± 80 


11.0 ±4.5 


-23 ± 12 


680 ± 80 


#3a 


7.5 ± 3.6 


-26 ± 1 1 


680 ± 80 


7.5 ±3.7 


-26 ± 12 


680 ± 80 


7.6 ±3.7 


-26 ± 12 


680 ± 80 


#3b 


16.6 ±8.6 


-95 ± 12 


690 ± 30 


17.0 ± 8.9 


-95 ± 12 


690 ± 30 


15.8 ±8.5 


-93 ± 12 


730 ± 20 



Notes, dp is the distance 



p the time, and v p the relative velocity at periastron of HE 0437-5439 with respect to the center of the LMC (runs #la, 
#2, and #3a) or to the GC (runs #lb and #3b). Numbers are mean value ± standard deviation cr. The initial conditions of runs #1 are those of 
Brown et al. (2010). In runs #2 and #3a, Kroupa & Bastian s ( 1997) LMC proper motions are used. In runs #3, proper motions of HE 0437-5439 
are increased by their systematic errors (CTI correction, see text for details). 
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Fig. 8. Histograms showing the distribution of distances d p , times T p , and relative velocities v p at periastron of HE 0437-5439 with 
respect to the LMC (upper panel, run #3a) and the GC (lower panel, run #3b) for Model III. Mean values and standard deviations are 
given in Tablep The gray-shaded areas mark regions with d p < R^ c = 3 kpc, d p < R^ c - lOkpc, and 15 Myr < |T p | < 21 Myr. 
The last is the lifetime of HE 0437-5439 assuming a single-star nature (Przybi lTa et al.|20 08 ). 



Table 4. Kinematic input parameters. 



Parameter 




Value 






HE 0437-5439 


LMC 


a (J2000.0) 




4 h 38 m 12 s .8" 


5 h 2T?6 d 


8 (J2000.0) 




-54°33'12" D 


-69°52' rf 


distance (kpc) 




61 ±9* 


50.1 ± 2A e 


fi a cos 5 (mas yr^ 1 ) 


0.53 


± 0.25(stat) ± 0.33(sys) c 


2.03 ± 0.08' 


£i,s (mas yr -1 ) 


0.09 


±0.21 (stat) ±0.48(sys) f 


0.44 ± 0.05' 


VradCkms -1 ) 




723 ± 3" 


262.2 ± 3.4 rf 



Notes. Uncertainties are lcr. 



References. ia) 


Edelmann et al. (|2005k; (b) 


Przybilla et al. J2008); 


(c) Brown et al. 


( 


2010); van der Marel et al. (2002); le ' Freedman 


etal. (2001 J; [J> 


Kallivayalil et al. 1 2006 1. 



When hypervelocity stars were first discovered in 2005 
( Brown et al.|2005 i, dynamical ejection from the supermassive 
black hole at the GC ( |Hills|19 88 ) was supposed to be their only 
origin because these stars move so fast that they are unbound to 
the Galaxy. This scenario was challenged by the discovery of the 
third hypervelocity star HE 0437-5439 (Ed elmann et al.|2005) , a 
massive B star, because the travel time from the GC to its present 
position would exceed its lifetime by a factor of 4. Edelmann 
et al.| ( |2005| ) notice that the star is much closer to the Large 
Magellanic Cloud (LMC) than to the Galaxy and might therefore 
originate in the LMC. This was corroborated by a sophisticated 
differential abundance analysis of high-quality, high-resolution 
spectra by Przybilla et al. (2008 ), who shows that the abundance 
pattern is inconsistent with that of GC stars but consistent with 
that of LMC B stars to within error limits. A kinematic inves- 
tigation for the place of birth of HE 0437-5439 would require 
precise measurements of the star's, as well as LMC's, proper 
motion. While several studies have been carried out to mea- 
sure the proper motion of the LMC, no measurements have been 
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130 140 150 160 170 180 190 200 
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Fig. 5. Visualization of the correlation between and in 
Model II. The meaning of the cross, curves, and shaded regions 
is the same as in Fig. [3] except that the red dashed line defines 
loci of constant halo mass inside a sphere of radius 20 kpc. 




50 

«h (kpc) 

Fig. 7. Visualization of the correlation between Mj, and ah in 
Model III. The meaning of the cross, curves, and shaded regions 
is the same as in Fig. [3] except that the red dashed line defines 
loci of constant halo mass inside a sphere of radius 18.8 kpc. 



available for HE 0437-5439 until recently. Bro wn et al.] ( f20T0l > 
used two epochs of images taken with the Advanced Camera 
for Surveys (ACS) onboard HST to determine the proper mo- 
tion of HE 0437-5439 and argue that their measurement rules 
out a place of birth in the LMC at the 3<x level but is consis- 
tent with an origin from the GC. To remedy the time-of-fiight 
versus lifetime problem, HE 0437-5439 needs to be a blue strag- 



gler (Edelman n et al.|2005| l. This would imply that the progen- 
itor was a triple system from which a binary was ejected by the 
supermassive black hole. Later on its trajectory, the binary sys- 
tem merged to form the blue straggler (Perets 2009). However, 
the result of Brown et al. (2010 1 is strongly affected by the treat- 



ment of LMC motions, as well as the charge-transfer inefficiency 
(CTI) correction applied to the data. To demonstrate both issues, 
the kinematical analysis of HE 0437-5439 is revisited here con- 
sidering three different cases. 



1. 



To start with, trajectories for 10 4 Monte-Carlo realizations, 
which simultaneously account for statistical and systematic 
uncertainties in current positions and velocities via Gaussian 
distribution, are computed for HE 0437-5439 and LMC in- 
dependently using the input parameters listed in Table [4] 
The distributions of pericenter distance d v , its correspond- 
ing point in time T p , and relative velocity v p derived from 
the resulting sample of 10 8 combinations of orbit pairs are 
insensitive to the choice of the mass model (see run #la in 
Table B), and they confirm the results of Brown et al. (2010 1 
(dp = 13 kpc, T p = -23Myr) based on the Kenyon et al 



(2008) Galactic potential and identical input values. Since 
dp exceeds the radii of LMC's central region R^ c = 3 kpc 



as well as the outermost regions R^ c - 10 kpc ( |Brown et al. 
2010), an origin in the LMC seems unlikely in this context. 
Next, the LMC proper motions of Table [4] are repl aced by 
<ji a cos6,n 6 ) = ( +1.94±0.29,-0.14±0.36)masyr 1 ( |Kroupa 
& Bastian]|1997 1 to explore their influence on the outcome. 
Performing the same Monte-Carlo method as before, all 
three models give smaller pericenter distances (see run #2 
in Table[5]) and thus closer encounters of the two objects. 
Assuming the CTI of ACS to be rising linearly in time, 
|Brown et al.| ( |2010| applied 55% of the epoch-2 CTI correc- 
tion to their epoch-1 images. According to Massey (2010), 
however, there was a dramatic increase in the CTI between 
the two epochs that lead to an overcorrection in the epoch-1 
data by Brown et al. implying proper motions of HE 0437- 
5439 that were larger than stated in Brown et al. (2010[). The 



corresponding effects are roughly estimated here by adding 
the systematic errors to the mean value and omitting them af- 
terwards, i.e., by using (//„ cos 6, us) - (+0.86±0.25, +0.57 ± 
0.21) mas yr _1 as input for the Monte-Carlo calculation. The 
resulting distributions (see run #3a in Table |5]) are still al- 
most model independent and visualized for Model III in the 
upper panel of Fig. [8] In all three models, about 10% of all 
orbit pairs yield pericenter passages within the central re- 
gion, i.e., dp < ^cen C ' or ^6% in the outermost regions of 
the LMC (dp < R^ c ). Out of these, 2% (21 %) have shorter 



flight times than the star's age of 18 ± 3 Myr (Przybilla et al 
2008). Thus, 0.2% (16%) of the trajectories are consistent 



with an origin in the inner (outer) LMC without invoking 
additional requirements, such as a blue straggler nature. The 
decrease in pericenter distances due to CTI effects (Table [5] 
run #2 versus run #3a) is stronger than when they are due to 
a change in LMC proper motions (Table [5] run #la versus 
run #2). 

For comparison, we now consider the GC as the place of ori- 
gin. Using the input values of Table [5] HE 0437-5439 passed 
95 + 14 Myr ago within d p = 14.8 ± 10.1 kpc of the GC (see 
Model III of run #lb in Table |5). While the travel time is in 
good agreement with Brown et al. (2010 1 (r p = -98 Myr), it 
is difficult to compare d p since it is not given in Brown et al. 
( |2010| ). Nevertheless, the upper lefthand panel of Fig. 3 in Brown 
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|et al.| ( 20T0| suggests that substantially more than 68% of all tra- 
jectories have Galactic plane-crossing locations below 15kpc. 
Because the Galactic plane-crossing location is an upper limit 
for a!p, this fraction increases further when we consider the peri- 
center distance. However, this contradicts our result that less 
than 60% of all orbits pass the GC within 15 kpc. As seen from 
run #3b in Table [5] the situation does not change significantly 
when accounting for the CTI overcorrection. The respective dis- 
tributions from the 10 4 trajectories of HE 0437-5439 are shown 
in the lower panel of Fig. [8] whereby only ~2% (27%) of the tra- 
jectories have pericenter passages within 3 kpc (10 kpc). All of 
them exceed the stellar lifetime and therefore require HE 0437- 
5439 to be a blue straggler. 

As a consequence, ruling out an origin in the LMC in fa- 
vor of the GC is disputable and the question of HE 0437-5439's 
birthplace remains open. 

7. Conclusions 

The motions of stars provide important information about the 
mass distribution in various components of the Galaxy. In par- 
ticular, they trace the dark matter. Oort s (1932) measurements 
of stellar motions were the first to hint at the presence of dark 
matter in the Milky Way. Soon, the Gaia mission will provide 
velocity information of a huge number of stars and satellite 
galaxies with unprecedented precision. Analytical mass mod- 
els of the Galaxy are utilized to calculate the orbits of stars. 



Allen & Santillan ( 1991 ) derived such a model by making use 
of observational constraints such as the Galactic rotation curve, 
the distance to the GC, and the local mass density. During the 
past twenty years, observational data have greatly improved. 
Therefore, it was time to revisit the Galactic gravitational po- 
tential of Allen & Santillan ( 1991 ) and update its parameters by 
making use of the latest observations. The halo mass is hereby 
constrained by the most extreme halo star discovered so far 



Therefore, we calculated trajectories for the star by making use 
of the three mass models under study. Using the same input data 
as Brown et al.| ( |2010"| , we could confirm their results regarding 
the LMC, while those for the GC show discrepancies that are in- 
dependent of the applied mass model. Our trajectories pass the 
GC within a considerably larger distance, showing that an ori- 
gin in the GC is much less likely than suggested by |Brown et al.| 
(2010). Moreover, we investigated the influence of LMC proper 
motions and inspected systematic errors in the star's proper mo- 
tion measurements stemming from the CTI overcorrection in the 
ACS images. The latter have not been considered by |Brown et al. 
(2010) but turned out to be crucial. Our calculations show that 
an origin of HE 0437-5439 in the LMC is at least as likely or 
unlikely as an origin in the GC. Presently available proper mo- 
tions are therefore inconclusive, and improved measurements are 
mandatory for settling this issue. 
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Appendix A: Equations of motion 

The motion of test particles in an axisymmetric gravita- 
tional potential 0(r, z) can be conveniently derived from the 
Lagrangian X. of the system: 

£(r,z,r,j>,z) = ^(r 2 + (r0) 2 +r)-(D(r,z). (A.l) 



( Przybilla et al 



2008}. For comparison., two other widely used introducing canonical momenta p r = 8£/dr = r, P4> = d£/dj> = 

r 2 (j), and p z = d£,/dz = z, the corresponding Lagrange's equa- 



halo mass distributions - the truncated, flat rotation curve model 
( |Wilkinson & E vans 1999} and a model derived from numeri- 
cal cosmological simulations (Navarro et al. 1997 1 - are fitted as 
well. 

All three Milky Way mass models are analytical, simple, and 
equally capable of reproducing their imposed restrictions. Major 
discrepancies only become apparent at Galactocentric distances 
greater than ~ 15 kpc, which is the region where observational 
constraints are rare and the halo component dominates. In par- 
ticular, the depth of the gravitational potential is very sensitive to 
the form of the dark matter halo seriously affecting for instance 
predictions of the local escape velocity. The significantly differ- 
ent shapes of the halos allow systematic uncertainties in kine- 
matic investigations to be estimated by comparing the results of 
orbit computations performed in the three models separately. 

Finally, the enigmatic hypervelocity star HE 0437-5439 is 
re-investigated by inspecting its trajectory in the context of the 
updated potentials. There are several pros and cons of whether 
the object originated in the LMC as proposed by the discoverers 
of the star (Edelma nn et al.|20 05 ) or from the GC, which is the 
suspected place of origin of all hypervelocity stars ( |Brown et al. 
|2005| ). The problem with the latter is that the travel time to the 
GC exceeds the stellar lifetime and that the chemical composi- 
tion of HE 0437-5439 differs from what is known about stars in 
the GC. A birthplace in the LMC, on the other hand, would be 
consistent with both. However, based on their own measurement 



tions - written as a system of first-order differential equations 
read as 



r = p r 



Pi 

2 ' 



r 

i = Pz 



P4> 

Pr = — 

*0 



(A.2) 
(A.3) 
(A.4) 



GM b r 



GM A r 



(r 2 + z 2 + 



b h 2 ) 



3/2 



+ (a d + Vz 2 + ^d 2 ) 



2\ 3 /2 



dDhlR) r 
dR Vr 2 + z 2 



P4, = , 



(A.5) 



(A.6) 



GM b z 



GM d z 



+ 1 



3/2 



of the star's proper motion, Brown et al. (2010) claim that the 
LMC - in contrast to the GC - is ruled out as the place of origin. 



(r 2 +z 2 + b h 2 ) 

d$hCR) z 
dR Vr 2 + z 2 



2\ 3 / 2 



(A.7) 
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whereby 

ggh(g) 
dR 



r-i 



GM h (f) 
ahR 1 + (Rj- 1 

GM h itT 



R 2 



GM h 



ifR< A & Model I, (A. 8a) 



ifR > A & Model I, (A.8b) 



y/R 2 + a h 2 R 
GM h I R 
^ ln,1 + 



GM h 



a h ) ( a h + R)R 



if Model II, 



if Model III. 



(A.8c) 



(A.8d) 



Expressing lengths in kpc, masses in galactic mass units M ga i, 
and time in Myr, the gravitational constant G has to be replaced 
by the factor ~ 1.04598 x 1(T 4 thereby converting units from 
100 km 2 s~ 2 to kpc 2 Myr" 2 . 
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